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ABSTRACT 

We investigate empirical scaling relations between the thermal Sunyaev-Zeldovich effect (SZE) and 
cluster mass in simulated clusters of galaxies. The simulated clusters have been compiled from four 
different samples that differ only in their assumed baryonic physics. We show that the strength of the 
thermal SZE integrated over a significant fraction of the virialized region of the clusters is relatively 
insensitive to the detailed heating and cooling processes in the cores of clusters by demonstrating that 
the derived scaling relations are nearly identical between the four cluster samples considered. For our 
synthetic images, the central Comptonization parameter shows significant boosting during transient 
merging events, but the integrated SZE appears to be relatively insensitive to these events. Most 
importantly, the integrated SZE closely tracks the underlying cluster mass. Observations through 
the thermal SZE allow a strikingly accurate mass estimation from relatively simple measurements 
that do not require either parametric modeling or geometric deprojection and thus avoid assumptions 
regarding the physics of the ICM or the symmetry of the cluster. This result offers significant promise 
for precision cosmology using clusters of galaxies. 

Subject headings: galaxies:clusters:general-cosmology:observations-hydrodynamics- 

methods: numerical-cosmology: cosmic microwave background 



1. INTRODUCTION 

Observations of galaxy clusters are increasingly used 
to constrain cosmological parameters. Reliable determi- 
nations of the cluster mass, abundance of clusters as a 



function of re dshift (e.g. iHaiman et all (j2001I) ). and gas 
fraction (e.g. Vikhli nin et alJ (|2003l )) lead to indepen- 
dent constraints on Qm, ^6? an d th e dark ener g y equa- 
tion o f state param e ter w (see e.g. iPenl (|1997l) . Sasaki 
(1996j), lAllen et all ((2004)). In most cases, extraction 
of cosmological parameters from clusters requires a high 
precision measurement of cluster mass, often for a large 
sample of clusters. It is therefore critical to understand 
whether existing or planned X-ray or SZE observations 
can yield sufficiently precise measurements to permit ac- 
curate determinations of cosmological parameters. 

Historically, cluster gas mass is often deduced from the 
X-ray temperature using an empirical relation between 
mass and spectral temperature with relatively small sca t- 
ter for high mass clusters (e.g. iSanderson et alJ (|2003l) ). 
The X-ray surface brightness of a cluster results from the 
emission along the line of sight, 



In contrast, the Sunyaev-Zeldovich effect (SZE) arises 
from the inverse Compton scattering of cosmic mi- 
crowave background photons by electrons in the hot 
plasm a of clusters of galaxies ((Sunvaev fc Zeldovichl 
1972). The strength of the thermal SZE is proportional 
to the Compton parameter, y, which for non-relativistic 
electrons is essentially the integral of the gas pressure 
through the cluster 
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where the proportionality holds if A includes only 
bremsstrahlung emission; which is a reasonable approx- 
imation for massive clusters of galaxies. In reality, the 
emission function will include both resonance lines and 
free-free emission. 
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For ver y hot clusters, re lativist ic corrections must be in- 
cluded ([Itoh et all 19981) . The central value of the Comp- 
ton y parameter will be referred to as yo. We also con- 
sider integrals of the Compton parameter over a finite 
projected radius; for example 2/500 is the integral of y 
over a disk with of radius rsoo corresponding to an over- 
density of 500 relative to the critical density. 

From the above discussion, these two cluster observ- 
ables depend differently on the state of the cluster gas. 
The X-ray emission is more sensitive to the density, and 
so one might expect X-ray observations to provide a bet- 
ter constraint on the cluster mass. However, this strong 
density dependence means that the emission is domi- 
nated by the cluster core, where the gas is subject to 
many complicating effects, including both cooling and 
heating. 

The two methods also have very different selection 
biase s. By its nature, th e SZE is redshift indepen- 
dent ([Carlstrom et al.ir2002l) and observations of clusters 
through this window hold great promise to test cosmo- 
logical models. Since the thermal SZE measures the in- 
tegral of the gas pressure and assuming that the gas out- 
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side the core is in approximate hydrostatic equilibrium 
within the dark matter potential, one expects that SZE 
observations should provide robust measures of the clus- 
ter mass - independent of the complicated interplay of 
physical mechanisms that regulate t he thermal state of 
the cluster gas (e.g. Bar bosa et alJ (|l996DlBenson et all 
(2004). 

To investigate this expectation, we have constructed 
four catalogs of simulated clusters of galaxies that dif- 
fer only in the assumed input physics for the baryonic 
component. Each sample contains ~ 100 clusters at 
the present epoch in the mass range from 10 14 M 
to 2 x 10 15 Mq and ~ 10 clusters more massive than 
10 14 Mq at a redshift of 2, roughly corresponding to the 
expected sensitivity of upcoming SZE telescopes. 

In this Letter, we use ideal synthetic observations of 
the thermal SZE and X-rays to measure how cluster ob- 
servables scale with cluster mass. Since no instrumental 
effects are included, this study indicates the minimum er- 
rors one can expect in these relations. Evaluating typical 
scaling relations for clusters reveals that the integrated 
SZE is a surprisingly accurate and simple method for de- 
termining the mass of clusters. Thus, SZE observations 
have a high potential for precision measures of funda- 
mental cosmological parameters. 

2. NUMERICAL SIMULATIONS 

Our simulations use the hybrid Eulerian adapt iv e mes h 
refineme nt /N-body code En zo (|Q'Shea et al.1 ((2004); 
http://cosmos.ucsd.edu/enzo) to evolve both the dark 
matter and baryonic fluid in the clusters, utilizing the 
piecewise parabolic method (PPM) for the hydrodynam- 
ics. With up to seven levels of dynamic refinement in 
high density regions, we attain spatial resolution up to 
~ 16 h~ x kpc in the clusters. We assume a concor- 
dance ACDM cosmological model with the following pa- 
rameters: Q h 0.026, Q m = 0.3, Q A 0.7, h 0.7, 
and erg = 0.928. Refinem ent of high density regions is 
performed as described in iMotl et all (|2004ft . We have 
constructed a catalog of AMR refined clust ers identified 
in the simulation volume as described in iLoken et alJ 
(2002). The four catalogs progressively include addi- 
tional physical processes in the calculation. Our baseline 
catalog was run with only adiabatic physics, a second 
includes the effects of radiative cooling, a third adds a 
model for the loss of low entropy gas to stars, and a final 
simulation adds a moderate amount of supernova feed- 
back to the star formation prescription. 

The loss of energy to radiation is calculated from 
a tabulated cooling curve deri ved from a Ravmond- 
Smith plasma emission model (Brickhous e et alJ 11995) 
assuming a constant metallicity of 0.3 relative to so- 
lar. For the simulation runs including star formation, we 
have incorporat ed the star formation prescription from 
ICen fc OstxikeTl fiflfll. 

Briefly, for every most refined cell, if the gas is con- 
tracting, cooling rapidly, and contains more than a Jeans 
mass of material, a new star particle is created with a 
mass raft rj At/td yn and this amount of mass is removed 
from the fluid. The star formation rate is thus coupled 
to the local dynamical time, td yn , while r] parameterizes 
the efficiency of star formation and At is the simulation 
timestep increment. In the simulation including feedback 
from stars, the new star particle begins to deposit energy 
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Fig. 1. — Images of the projected, emission- weighted temperature 
map (left column), normalized X-ray surface brightness (middle 
column) and the Compton parameter, y (right column). Each row 
displays the same cluster evolved in the indicated physical limit 
at the present epoch. The field of view is 5 Mpc on a side. 
The color table for each image is set by the range of values in 
the adiabatic realization. This particular cluster has a total mass, 
M2oo(~ M viria i) ~ 2 x 10 15 Mq. In reality, clusters should be 
intermediate between the limiting cases of radiative cooling only 
(which produces too many clusters with cool cores) and the star 
formation feedback model (which produces too few). 

in the fluid to simulate the explosion of prompt, type II 
super novae. 

3. RESULTS 

Figure 1 depicts typical simulation results from the 
four cluster samples that we have considered. The most 
dramatic difference between the samples is the rich ar- 
ray of cool substructures in the cooling only simulation. 
These cool cores and their associated shock features are 
very conspicuous in the temperature map and X-ray im- 
age. While there is significant variation between clusters 
when viewed through their X-ray emission directly or 
through the derived temperature maps, the thermal SZE 
yields a quite similar appearance for the clusters - even 
in the extreme case of cooling only. 

3.1. The Effect of Cluster Mergers 

The bias introduced by structure formation itself raises 
potential concerns, independent of the complexities of 



3 



TABLE 1 

Scaling Exponent for ysoo - 



M Relation, z=0 



TABLE 2 

Accuracy of Mass Estimation, z=0 



Simulation 


a 




Method 


Median M est /M t rue 


+ 1<7 


-la 


+80% 


-80% 


Adiabatic 


1.59 


0.021 


2/500 - M 


0.97 


1.00 


0.93 


1.13 


0.86 


Radiative Cooling 


1.71 


0.031 


T x -M 


1.04 


1.11 


0.92 


1.33 


0.74 


Star Formation 


1.60 


0.027 


^X,500 — M 


0.87 


1.03 


0.75 


1.46 


0.62 


Star Formation with Feedback 


1.61 


0.024 


yo-M 


0.96 


1.14 


0.82 


1.53 


0.67 



the detailed physics governing energy balance in the clus- 
ter gas. Generally, we find that the central value of the 
Compton parameter (yo) can be boosted by up to a fac- 
tor of ~ 20 in head on collisions between equal mass sub- 
clusters. These major mergers can drive the evolution of 
the central SZE signal over time periods of 1 — 2 td yn or 
~ 2 Gyr. Furthermore, even minor mergers (which are 
more common) with a mass ratio of ~ 0.1 can produce 
significant boosts in yo by a factor of ~ 2. 

On a cluster by cluster basis, merger boosts will bias 
measurements derived from the cluster properties to, 
for example, overestimate the mass of the cluster cor- 
responding to a given SZE signal. In a survey of clus- 
ters, merger boosts will scatter systems with equilibrium 
values below the detection threshold (which are more nu- 
merous) into the sample. 

While mergers significantly impact yo, they do not im- 
pact the value of 2/500 as drastically. Measuring the in- 
tegrated SZE out to a significant radius has the effect 
of de-emphasizing the core properties of the cluster, and 
so the integrated SZE is less sensitive to the dynamical 
state of the cluster, as well as to the details of the cluster 
physics as we show below. 

3.2. Thermal SZE Scaling Relations 

The comparison of the 2/500 — ^500 relation with the 
yo — M relation is shown in the upper panels of Figure 
2 for the star formation with feedback sample. The yo 
scaling shows a much larger scatter than 2/500, and so 
mass estimation using the central value of the Compton 
parameter is correspondingly poorer. To illustrate the 
reason for this discrepancy, we show in the bottom pan- 
els of Figure 2 both the central pressure in each cluster, 
and the integrated pressure inside a sphere of radius r$oo 
as functions of mass, yo is essentially a measure of the 
pressure integrated along the line of sight through the 
cluster center. However, the integrated SZE, 2/500, mea- 
sures the projected integrated pressure inside r^oo- It 
is clear from these plots that while the central pressure 
is weakly correlated with cluster mass, the value of the 
pressure integrated out to a large radius is a very good 
predictor of the mass. This difference is a result of a vari- 
ety of effects, including both mergers and non-adiabatic 
physics, which dominate the cluster cores. Observing the 
integrated SZE signal more effectively samples the clus- 
ter potential depth by measuring the integrated pressure 
out to a large radius, leading to a tight correlation of the 
strength of the signal with cluster mass. Looking back in 
time through our catalogs of simulated clusters, we find 
that all four samples produce similar predictions for the 
scaling between 2/500 and M500. These predictions are 
roughly consistent with the expectations of self-similar 
scaling for clusters, up to a redshift of 1.5 for which we 
have a sufficient number of clusters in our simulations. 
Table 1 shows the scaling exponent and the variance for 



each of the four samples at z = 0, where the scaling 
relation is written as 
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The values of the exponent a s hown in the table ar e 
consistent with previous work bv Ida Silva et all ((2004) . 
While the detailed simulation of clusters incorporating 
all relevant physics remains an unfulfilled challenge, it 
is reassuring that all models considered by us to date 
predict similar properties for the integrated SZE. This 
results from the weak dependence of the integrated SZE 
on the state of the cluster core. Therefore, the inter- 
pretation and modeling of SZE observations of clusters 
should be relatively immune to theoretical uncertainties 
and SZE surveys will provide excellent opportunities for 
constraining cosmological parameters. 

3.3. Cluster Mass Estimates 

The relationship between X-ray spectral temperature 
and cluster mass (M — Tx) is well documented and a 
similar relation is noted in cosmological simulations (e.g., 
iBrvan Normanl (|1998t )V This relationship is often used 
to infer cluster masses. For comparison to the 2/500 — M 
relation studied here, we have measured several other 
scaling relations, including M — Tx, for the simulated 
clusters. Tx is in this case calculated as the X-ray emis- 
sion weighted average temperature inside a projected ra- 
dius of r5oo- We then use the best-fit scaling relations 
derived from the star formation with feedback simula- 
tion to determine the total cluster mass. Table 2 shows 
the ratio of estimated mass to true mass for four different 
scaling relations at z=0. 

We find that the integrated thermal SZE signal within 
a characteristic radius of r$oo is the best estimator of 
the true cluster mass among these scaling relations. For 
example, 80% of clusters have estimated masses that lie 
within +15% to —10% of the true cluster mass and this 
accuracy is nearly constant with redshift back to at least 
z = 1.5. The 2/0 value is a much poorer estimate of the 
true cluster mass, as is expected from the relative scatter 
of yo — M shown in Figure 2. Core effects also have a 
strong impact on the X-ray luminosity, because of its 
enhanced density dependence. Therefore the scaling of 
X-ray luminosity with cluster mass is the poorest of these 
scaling relations. 

The M — Tx relation for clusters in the simulation gen- 
erate s relative errors of ±30% in mass for 80% of clus- 
ters. llRasia et"aTI (|2005l) show that, from simulations, a 
la scatter of 30% in the M — Tx relation is reduced to 
16% when using their spectroscopic-like temperature in 
place of the emission weighted temperature. A similar 
drop in the scatter of the mass estimates quoted here 
for 80% of the clusters would make the M — Tx rela- 
tion roughly equivalent to 2/500 — M in predictive power. 
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Fig. 2. — Upper: The scaling relations between projected yo and 2/500 and the total cluster mass within the same radius at the present 
epoch for the star formation with feedback cluster sample. Two randomly chosen, orthogonal projections for each cluster are plotted as 
individual points and the catalog contains ~ 100 clusters at this epoch in the mass range 10 l4 Mq < M200 < 2 x 10 15 Mq. The best 
fit relations are plotted as solid lines. Lower: Central pressure and pressure integrated inside sphere of radius rsoo plotted against cluster 
total mass. 



However, this analysis includes no instrumental effects; 
a detailed study of the additional error introduced using 
real observations remains to be done in order to gauge 
the effectiveness of these methods. Certainly, measuring 
the value of 7/500 requires no assumptions about the state 
of the cluster gas, nor any model fitting, whereas X-ray 
temperature is generally fit to a model assuming uniform 
temperature and metallicity within the gas. SZE obser- 
vations also can be done for a larger number of clusters, 
and at higher redshifts, than X-ray observations. In any 
case, of the methods tested here, the integrated SZE scal- 
ing relation is superior at estimating cluster masses. 

4. CONCLUSIONS 

The integrated SZE in clusters of galaxies is not 
strongly dependent on the details of cluster physics or 
variation in the cluster core in contrast to the X-ray lu- 
minosity. 

While the value of yo experiences strong "boosting" 
due to mergers, the value of 7/500 appears to be much less 



sensitive to transient events. Measuring the SZE out to 
larger radii has the effect of smoothing out the impact of 
mergers. 

The 7/500 — M relation is the simplest, least model de- 
pendent, and most accurate measure of cluster mass. 
An unprecedented number of clusters will be observed 
through the SZE in the near future. Clusters observed 
through the SZE therefore show great promise as preci- 
sion cosmological probes. 
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